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IS IT NECESSARY TO GO BEYOND THE PONCTUAL MASS APPROXIMATION 
FOR TIDAL PERTURBER IN CLOSE SYSTEMS? 

Mathis, S.^ and Le Poncin-Lafitte, C.^ 

Abstract. With the discoveries of very close star-planet systems with planet orbiting sometimes at several 
star radius but also with well-known situations in our solar system where natural satellites are very close to 
their parent planet the validity of the ponctual mass approximation for the tidal perturber (respectively the 
parent star or planet when we study the close planet or natural satellite dynamics) has to be examined. In 
this short paper, we consider this problematic using results coming from a complete formalism that allows to 
treat the tidal interaction between extended bodies. We focus on its application to a simplified configuration. 

1 Context 

In celestial mechanics, one of the main approximation done in the modelling of tidal effects (star-star, star- 
innermost planet or planet-natural satellites interactions) is to consider the tidal perturber as a point mass 
body. However a large number of extrasolar Jupiter-like planets orbiting their parent stars at a distance lower 
than 0.1 AU have been discovered during the past decade (Mayor et ai, 2005). Moreover, in Solar System, 
Phobos around Mars and the inner natural satellites of Jupiter, Saturn, Uranus and Neptune are very close to 
their parent planets. In such cases, the ratio of the perturber mean radius to the distance between the center 
of mass of the bodies can be not any more negligeable compared to 1. Furthermore, it can be also the case 
for very close but separated binary stars. In that situation, neglecting the extended character of the perturber 
have to be relaxed, so the tidal interaction between two extended bodies must be solved in a self-consistent 
way with taking into account the full gravitational potential of the extended perturber, generally expressed 
with some mass multipolc moments, and then to consider their interaction with the tidally perturbed body. In 
the litterature, not so many studies have been done (Borderies, 1978-1980; Ilk, 1983; Borderies & Yoder, 1990; 
Hartmann, Soffcl & Kioustelidis, 1994; Maciejewski, 1995). 



2 General formalism 

Several years ago, Hartmann, Soffel & Kioustelidis (1994) introduced in Celestial Mechanics an interesting 
tool, based on Cartesian Symmetric Trace Free (STF) tensors, to treat straighforwardly the couplings between 
the gravitational fields of extended bodies. These tensors are fully equivalent to usual spherical harmonics 
but in addition a set of STF tensors represents an irreductible basis of the rotation group SOS (Courant & 
Hilbert 1953). It means that using algebraic properties of STF tensors, these objects become a powerful tool to 
determine the coupling between spherical harmonics in an elegant and compact way. However, as these tensors 
are not widely used in celestial mechanics, Mathis & Le Poncin-Lafitte (2008) (hereafter MLP08) first recall 
their definition and fundamental properties and stress their relation with usual spherical harmonics. Then, the 
multipole expansion of gravitational-type fields is treated. First, the well-known external field of such body is 
derived using STF tensors; classical identities are provided. Next, the mutual gravitational interaction between 
two extended bodies and the associated tidal interaction are derived. We show how the use of STF tensors leads 
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to an analytical and compact treatment of the coupling of their gravitational fields. We deduce the general 
expressions of tidal and mutual interaction potentials expanded in spherical harmonics. Using the classical 
Kaula's transform (Kaula, 1962), we express them as a function of the Keplerian orbital elements of the body 
considered as the tidal perturber. These results are then used to derive the external gravitational potential 
of such tidally perturbed extended body. Introducing a third body, its mutual interaction potential with the 
previous tidally perturbed extended body is defined that allows us to derive the disturbing function using the 
results obtained with STF tensors and the Kaula's transform. The dynamical equations ruling the evolution of 
this system are obtained. 

3 Case of an extended axisymmetric deformed perturber 

Here, our goal is to quantify the term(s) of the disturbing function due to the non-ponctual behaviour of the 
perturber B and to compare it to the one in the ponctual mass case. To achieve this aim, some assumptions 




Fig. 1. Classical tidal dynamical system. The extended body B is tidally disturbing the extended body A which adjusts 
itself with a phase lag St^ due to its internal friction processes. The dynamics of B is then studied. 57a and ub are 
respectively the spin frequency of A and the mean motion of B. 

are done. First, we adopt the quadrupolar approximation for the response of A to the tidal excitation by B; 
thus, using the calculations of MLP08, the disturbing function is then given by 
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where Li = {toa; ^B; "Ib, <?}; indexes definition is given in MLP08. G is the universal constant of gravity. Mb 
is the mass of B and R\ is the radius of A. fc^ is the Love's number of A that gives the linear adiabatic response 
of A to the tidal perturbation exerted by B. Zta describes the dissipation of the tidal kinetic energy by viscous 
friction and thermal diffusion {v and K are respectively the (turbulent) viscosity and the thermal diffusivity) . 
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On the other hand, the gravific potential of B is expanded as (r) — G^;^>q Af^^ ^(b/tib {(^if) 
where the Mf^ are the gravitational multipole moments of B and (r, 9, ip) are the spherical coordinates which 
have the center of mass of B for origin, ae, ee and Ib are respectively the semi- major axis, the eccentricity and 
the inclination of the relative orbit of B. '^L.M,j,p,q is a function of the three other keplerian elements of this 
orbit, of the A angular velocity ^Ia and of its precession angle, ea is the obliquity of A. 7^^'^/^ and k^j are 
coupling coefficients, which are detailed in MLP08. d^M' ^Lj-.p ^^"^ ^L.p.q are respectively obliquity, inclination 
and eccentricity special functions. 



Since we are interested in the amplitude of "R-Li, we focus on its norm (|7?,ij|). On the other hand, as we 
know that the dissipative part of the tide is very small compared to the adiabatic one (cf. Zahn 1966), we can 
assume that I^TaI ~ 1 in this first step. 

Let us first derive the terms \'R-Li \ due to the non-ponctual terms of the gravific potential of B which have a 
non-zero average in time over an orbital period of B, (y^_-p)j,^ (r) = I/Tb /p ^ {t, r) dt that corresponds 

to the axisymmetric rotational and permanent tidal deformations (see Zahn 1977) (the same procedure can of 
course be applied to the non-stationnary and non-axisymmetric deformations, but we choose here to focus only 
on (V^_-p) to illustrate our purpose). Then, as the considered deformations of B are axisymmetric, we can 
expand them using the usual gravitational moments of B {Ji^ ) as 
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with p + M; — ^ j^iy where p and Af^ "^p are respectively the gravitational multipole moments 

induced by the permanent rotational and tidal deformations of B; Mb is the mass of B. Ml,m is the spherical 
harmonics (Il.a/) normalization constant. Then, we obtain 
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On the other hand, the term \TZli\ associated to Mb, namely the disturbing function in the case where B is 
assumed to be a ponctual mass, is given by: 
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focusing on the configuration of minimum 



energy. In this state, the spins of A and B are aligned with the orbital one so that = -^b = (that leads to 
J = TOA and p = (2 — mA + ^b) /2) and the orbit is circular (cb = 0). Then, we consider: 
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Using eqs. p.3ll3.4[) . we get its expression in function of Ji^ and of (/Jb/ib) 
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As it has been emphasized by Zahn (1966-1977), the main mode of the dissipative tide ruling the secular 
evolution of the system is tua = 2. We thus define £1-^ such that 
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which can be recast into 
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Finally, keeping only into account the quadrupolar deformation of B (J2), we get: 
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This gives us the order of magnitude of the terms due to the non-ponctual behaviour of B compared to the one 
obtained in the ponctual mass approximation. It is directly proportional to the squarred J2, thus increasing 
with (where en = ri|/J7^, fie being the angular velocity of B and Vic — \ ^ip^) the case of the rotation- 



induced deformation and with e\ (where £t — q{Rb/o.b) where q — AIj^/AI-q) in the tidal one, while it 
increases as (-Rb/qb)^- Therefore, as it is shown in Fig. the non-ponctual terms have to be taken into 
account for strongly deformed perturbers (J2 > 10"^) in very close systems {a^/R^ < 5) while they decrease 
rapidly otherwise. 
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